Observational constraints upon quintessence models arising from moduli fields 



O 

o 
o 

(N 

c 



(N 
> 

o 
o 
o 

Oh' 
I ' 

O 



13 



S.C. Cindy Ng 

Department of Physics and Mathematical Physics, University of Adelaide, 

Adelaide, S.A. 5005, Australia. 

(2 Jun, 2000; ADP-00-19/M92, [a.stro-ph/0004196; ) 

We study observational constraints on cosmological models with a quintessence arising from moduli 
fields. The scalar field potential is given by a double exponential potential V = Vo exp(— Ae^"*). 
After reviewing the properties of the solutions, from a dynamical systems phase space analysis, we 
consider the constraints on parameter values imposed by luminosity distances from the 60 Type IA 
supernovae published by Perlmutter et at, and also from gravitational lensing statistics of distant 
quasars. We also update the constraints on models with a single exponential potential V = Voe~ XK ^. 
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I. INTRODUCTION 

A cosmological constant A has been considered as the 
missing energy of the universe for a long time. Recently, a 
varying vacuum energy or "quintessence" [[lj has become 
popular as an alternative candidate. Since more parame- 
ters are involved in such models, they are more flexible in 
explaining some of the problems left by the cosmological 
constant model. One such problems is the "cosmic coin- 
cidence problem" Q : the missing energy and the matter 
energy densities decrease at different rates as the uni- 
verse expands, so it seems purely coincident that they 
are comparable today with the concordant values from 
several observational tests of fi^o ~ § and ^mo ~ \ |§ ■ 
Several candidates for quintessence fields have been pro- 
posed. Typically quintessence models possess attractor 
solutions with common evolutionary properties for a wide 
range of initial conditions. For example, in some cases 
the scalar field energy maintains constant ratio to the 
matter density at late time Q, where as in other models 
it tracks the dominant matter component in some other 
sense 

In this paper, I will study quintessences arising from 
string moduli. In string or Kaluza-Klein type models the 
moduli fields associated with the geometry of the extra 
dimensions may have effective potentials which depend 
exponentially on the moduli fields, due to the curvature 
of the internal spaces, or alternatively through the inter- 
action of the moduli with the form fields on the internal 
spaces (see, e.g., |?| and references therein). Single expo- 
nential potentials of the form 



V = Vne-^ 



(1) 



which give rise to a scaling solution have been well- 
studied in the literatures (see, e.g., and references 
therein). In this paper I will concentrate on double ex- 
pontial potentials of the form 



which typically arise from supersymmetry breaking via 
gaugino condensation |l(]-[l2|]. Furthermore, it has also 
recently been argued that if the superpotential of the 
moduli field in string theory is T-duality invariant, then it 
cannot be approximated by a single exponential function, 
but must de pen d on double exponentials pTsl. 

Binctruy [12], and later de la Macorra |jl3|], have ar- 
gued that scalar fields with double exponential poten- 
tials cannot act as quintessences for two reasons. Firstly, 
such models have a global attractor solution leading to a 
matter-dominated universe at late times. Secondly, near 
the attractor the equation of state is positive, w<j, > 
and approaching zero as t — > oo. This would appear to 
contradict the latest observational results from type IA 
supernovae that a quintessence is dominating over mat- 
ter and the universe is entering a phase of accelerated ex- 
pansion (UJlJ] , similar to inflation. The purpose of this 
paper is to show that for parameter values away from 
the attractor, there can exist models which are consis- 
tent with the observational tests, although some tuning 
is required for the universe not to have reached the at- 
tractor at present. This paper also updates the obser- 
vational constraints on the single exponential potential 
model, which have been given by Frieman and Waga [ fl6| . 

The organization of this paper is as followed: in sec- 
tion II, I will present a phase-space analysis of the dou- 
ble exponential potential model. In section III, I will 
discuss the numerical integration of the evolution equa- 
tions for both the single and double exponential potential 
models, and obtain values for Q^q and -ffo^o- I n section 
IV, I will constrain both models using the light-curve 
calibration luminosity distances of type IA supernovae 
(see [ pr|Jl7| and references therein) and the gravitational 
lensing statistics of high luminosity quasars by interven- 
and references therein). 



ing galaxies (see [18 



|19H and references therein). Similar 
constraints on other quintessence models have been pre- 
sented in Gafi§i. 



V = Vo cx-pi-Ae^t 



(2) 
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II. PHASE SPACE ANALYSIS FOR DOUBLE 
EXPONENTIAL POTENTIAL MODEL 

We begin by considering a universe which consists of a 
scalar held with the potential (||) and a barotropic fluid 
with equation of state P 7 = (7 — l)p 7 , < 7 < 2. For a 
spatially-flat Friedmann-Robertson- Walker (FEW) uni- 
verse, the governing equations are given by 



H 

Pi 



K ( 
-ZH{ Pl - 



P-y) , 



subject to the Friedmann constraint 



1 



V 



(3) 
(4) 
(5) 



(6) 



where an overdot denotes the ordinary differentiation 
with respect to the time t. In the above equations, 
k 2 = 87rG, and H = a/ a is the Hubble parameter. 

We follow the treatment of Copeland, Liddle and 
Wands M by defining the variables 



x = 



K0 



and introducing a third variable p2| ] 

In the case of the single exponential potential (|l|), \(</)) 
reduces to a constant. In terms of these variables the 
evolution equations (||) - @ become 

x' = -3x + Ay|y 2 + ^[2 a ; 2 +7(l- a; 2 -y 2 )] , (9) 
V' = -\\j\xy+ i -y[2x 2 + 1 {l-x 2 -y 2 )\ , (10) 



kVV 



(7) 



(8) 



A' = 2v / 3a;A 



(11) 



where a prime denotes a derivative with respect to the 
logarithm of the scale factor, N = ln(a). The system is 



bounded by the cylinder x 



y 



1. Assuming that A 



is positive, the system is confined to A > 0. The system 
is symmetric under the reflection y — > —y. 

As many properties of the solutions in the double 
exponential potential model can be related to the sin- 
gle exponential potential model, let us review the sin- 
gle exponential potential model as has been studied in 
the literature ( and references therein). Depend- 

ing on the values A and 7, there are up to five criti- 
cal points in the x - y plane, two of which are related 
to late times. For A 2 < 37, the late times attractor is 
a scalar field dominated solution = K 2 p^,/3H 2 — > 1 



with 70 ee (p^ + P$) I p0 — ► 4j- < 7. For A 2 > 37, the late 
times attractor is the scaling solution Q with 7^ = 7 and 
pj-, a constant. Fig. |l| shows the late times attrac- 
tors and how they evolve as A 2 increases from zero. As 
A 2 increases the attractor evolves from a stable node into 
a stable spiral at A 2 — 247 

the origin as A 2 
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FIG. 1. The evolution of the late times attractors of the 
single exponential potential model as A 2 increases. The scalar 
field dominated solution has x = A/%/6, y = yl — A 2 /6; the 
scaling solution has x = y // 3/2'j/\, y = -y/3(2 — 7)7/2A 2 . 
The scalar field dominated solution is a saddle point for 
A 2 > 37. 

In the double exponential potential model, we can 
identity critical points at finite A from the evolution 
equations (||) - (JlT|) . They are three discrete points at 
(x ci y ci \ c ) — (±1,0,0), (0,1,0), and a 1-parameter fam- 
ily at (0, 0, A c ) for A c arbitrary. The stability of the criti- 
cal points can be obtained by linearizing (g) - (|ll]) about 
the these points and solving for the eigenvalues of small 
perturbations. The critical points and the eigenvalues 
are listed in Table |]. 

TABLE I. The critical points at finite A and their eigen- 
values for the double exponential potential model. 



X c 


Vc 


A c 


Eigenvalues 


1 








3, 3, 2V3 


-1 








3, 3, -2y/3 





1 





-3, -3, 








< A c < 00 


-3/2, 3/2, 



The scalar field kinetic energy dominated solution 
(1,0,0) is an unstable node, while (—1,0,0) is a sad- 
dle point since trajectories are repelled in the x and y 
directions and are attracted in the A direction. The 
scalar field dominated solution (0, 1, 0) attracts a two- 
dimensional bunch of trajectories but is degenerate in 
the A direction. For the barotropic fluid dominated solu- 
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tion (0,0, A c ), < A c < oo, trajectories are attracted in 
the x direction, repelled in the y direction and degenerate 
in the A direction. 

To examine the critical points at infinite A, we use the 
Poincare sphere method by transforming 



A 



On the surface e = we find 



dx 



dy 
dr 



xy 



(12) 

(13) 
(14) 



where r is a new time coordinate defined by dr = XdN. 
Therefore, there are infinitely many critical points ly- 
ing on the e = surface with y c — and arbitrary 
x c € [—1,1]. The eigenvalues for the matrix of pertur- 



bations are 0, 



0. 



Therefore the critical points 

are degenerate in both the x and A direction. They repel 
or attract trajectories in the y direction depending on 
-1 < s c < or < it < 1 respectively. For x c = all 
the eigenvalues vanish, and higher-order perturbations 
show that it is a stable attractor. 

We will numerically integrate the evolution equations 
@ ~ (0) f° rwar ds and backwards from some initial 
point (xi,yi,\), in the case of a matter-dominated uni- 
verse 7=1. In Fig. U and we project the trajec- 
tories onto the x - y plane. When integrating back- 
words from the initial points (Xi,yi, Aj), the trajectories 
in Fig. H approach the critical point (1, 0, 0), and the tra- 
jectories in Fig. |^ approach the critical points (x c , 0, A c ), 
x c € [—1, 0), A — > oo. 

The scalar field is rolling down the potential at suffi- 
ciently late times so <fi and hence \(<p) is increasing with 
time. As A(<^>) increases with time, the behaviour of the 
solutions effectively approximates that of single exponen- 
tial potential models, but with the parameter A evolving 
in the direction shown in Fig. l|. For sufficiently small 
Xi the scalar field-dominated solution is dynamically im- 
portant and the trajectories get very close to the circular 
boundary. As A increases, the solutions approach scaling 
solutions (x = y for 7=1) before spiralling towards the 
origin. 



III. NUMERICAL INTEGRATION 

In order to determine il^o and Hgto, and the lumi- 
nosity distance - redshift relation in the next section, 
we ought to integrate the evolution equations (||) - (j|) 
numerically. In this section, we will discuss the numeri- 
cal integration for both the single and double exponen- 
tial potential models in a manner similar to a previous 
study on the model with a pseudo Nambu-Goldstone bo- 
son (PNGB) potential §l|- 
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FIG. 2. The projection of the space trajectories on the 
x - y plane for 7 = 1. The trajectories are for initial 
values (xi,yi,Xi): (0.1,0.1,0.3), (0.2,0.2,0.5), (0.3,0.3,1.0), 
(0.4, 0.4, 1.5). The dashed line is x = y. 
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FIG. 3. The projection of the space trajectories on the x 
- y plane for 7 = 1. The trajectories are for the initial val- 
ues (xi,yi,Xi): (-0.1,0.2,15), (-0.1,0.4,10), (-0.1,0.6,5), 
(-0.1, 0.8, 0.5). The dashed line is x = y. 



A. Single exponential potential 

Here we update the numerical integration result for the 
single exponential potential by Frieman and Waga |fl6| . 
We use an alternate set of dimensionless variables: 



u = — ,v = il„' l0 (l + z) ,w = Xk(, 
xio 



In 



Hi 



, (15) 



the variables u and v are the same as were used in [|2l|. 
In terms of these variables, the field equations become: 



-o—u + Xe 
-Ho 

H 



Ho 

w = Ait , 



(16) 

(17) 
(18) 



where a prime denotes a derivative with respect to the di- 
mensionless time parameter Hot. The Hubble parameter 



3 



is defined implicitly according to 



H ( 3 1 2 1 - 

Ho V 6 3 



(19) 



We shall begin the integration at last-scattering. We 
choose Hoti = 0. At this initial stage, H/Ho 3> 1 and 
the field evolution (|l6|) is overdamped by the expansion 
of the universe, driving u to zero. We therefore choose 



1/3. 



Ui = U. At last-scattering, 2;~11UU ana il^ ^l for 
f2m0~l) hence = 1100 is a reasonable value. Note 
that results of the integration do not change significantly 
if Hoti, Ui and i?j are altered to values within the same 
order of magnitude. The remaining parameters are Wi 
and A, which are determined by the physical origin of 
the model. 

The integration proceeds until a value H/Ho = 1 is 
reached. In Figs. || and [|, the contours of fl<po = 1 — v o 
and Hoto are displayed in the Wi - X parameter space. 
In view of recent estimates of the ages of globular clus- 
ters [|3|, a mean value of 12.8 Gyr for the age of the 
universe is arrive at. With ft, ~ 0.7 this would require 
Hoto ~ 0.9. A concordant value of CUq — 0.67 ± 0.05 
is given by several observational tests |H] . By comparing 
these two contour plots, we see that parameters in the 
region of small A are favoured. This region shall repre- 
sent models where the solutions approach the attractor, 
either the scalar field dominated solution (A 2 < 37) or 
the scaling solution (A 2 > 37), only recently. Note that 
a seperate bound of A 2 > 20 has been obtained from nu- 
cleosynthesis by assuming that the late times attractor 
is the scaling solution, and the attractor is approched 
within a few expansion times of the end of inflation . 
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FIG. 5. Contours of constant Hoto in the Wi — A plane for 
the single exponential potential model. 



qo < corresponds to a universe whose expansion is ac- 
celerating at the present epoch. 



B. Double exponential potential 

With the double exponential potential model, we again 
use the variables u and v as defined in (|l^) , and redefine 
w as 

w = V2K(p + hxA . (21) 
In terms of these variables the field equations become 
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FIG. 4. Contours of constant fl^o in the Wi - X plane for 
the single exponential potential model. 

Another useful quantity which can be determined from 
the numerical integration is the deceleration parameter 



3 1 



(20) 



-3^-u + V2Vie w exp(-e u ') 
Ho 

H 



-no 



V2u , 



(22) 

(23) 
(24) 



where we have replaced Vo with a dimcnsionlcss param- 
eter V\ — t^j^r, and the Hubble parameter is 

£1 n 



H 
H~o 



l 7 u 2 + \v x exp (-e w ) 
6 3 



(25) 



We choose the same HoU, Ui, and Vi as above. The pa- 
rameters Wi and V\ determine the phyisical origin of the 
model. The contours of fl^o an d Hnto are displayed in 
the Wi - Vi parameter space in Fig. 0and|^. By compar- 
ing these two contour plots, the interesting region would 
seem to be the bottom left-hand corner of the parameter 
space. The interesting region corresponds to the scalar 
fields being still nearly frozen to their initial states, or 
having become dynamical and starting to roll down the 
potential slope only recently. 
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FIG. 6. Contours of constant Q^o in the Wi - Vi plane for 
the double exponential potential model. 




FIG. 7. Contours of constant Hoto in the u>i - Vi plane 
for the double exponential potential model. 



IV. OBSERVATIONAL CONSTRAINTS 

A constraint on the single exponential potential model 
has been obtained by Frieman and Waga using five 
type IA supernovae (Sne IA). In this section, we will con- 
strain both the single and double exponential potential 
models using a larger set of Sne IA, together with a sep- 
arate constaint using the gravitational lensing statistics 
of high luminosity quasars. 

Empirical calibration of the light curve - luminosity re- 
lationship of Sne IA provides absolute magnitudes that 
can be used as distance indicators. The Supernova Cos- 
mology Project leaded by Perlmutter et al. and the 
High Redshift Supernovae Search Team leaded by Riess 
et al. [|l5| are two different groups have been working ex- 
tensively in searching for high redshift Sne IA. I will use 



the larger available data set of the two, namely the 60 
Sne IA published by Perlmutter et al. 10 . The lumi- 
nosity distance - redshift relation so obtained provides 
a good observational constraint on quintessence models 
provided that there is no intrinsic evolution of the peak 
luminosities of the Sn I A sources. The issue of how the 
constraints are altered in the presence of such evolution 
has been discussed in fi~7|j2l"|| . 

Gravitational lensing of distant quasars due to galax- 
ies along the line of sight provides another relatively sen- 
sitive constraint Jjl| on the quintessence models. The 
statistics of abundances of multiply imaged quasars and 
observed separations of the images to the source can be 
used to estimate the distances to the quasars. We follow 
the calculation as described by Waga and Miceli |n| , who 
used a total of 862 (z > 1) high luminosity quasars plus 
5 lenses from seven major optical survey. 

We display the results of both tests as 68.3% and 95.4% 
joint credible regions of the parameter spaces described 
in the previous section. Figs. ^ and ^ display results for 
the parameter spaces of the single and double exponen- 
tial potential models respectively. For both models, the 
interesting regions that give rise to expected values of 
^00 ~ 0.67 and H t ~ 0.9 are well within the 68.3% con- 
fident region of both the Sne IA and the gravitational 
lensing statistics tests. The regions to the left of the 
qo = contours correspond to models that give rise to 
accelerated expansion. 
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FIG. 8. The 68.3% and 95.4% joint credible regions on 
the Wi - X parameter space of the single exponential potential 
model, for the Sne IA test (dashed) and the gravitational 
lensing statistics test (dotted). Over plots are the contours of 
Hoto = 0.9, Qtpo = 0.67, and qo = 0. The region allowed at 
95.4% by both observational tests is shaded. 
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FIG. 9. The 68.3% and 95.4% joint credible regions on the 
Wi - Vi parameter space of the double exponential potential 
model, for the Sne IA test (dashed) and the gravitational 
lensing statistics test (dotted). Over plots are the contours of 
Hoto = 0.9, fi</,o = 0.67, and qo = 0. The region allowed at 
95.4% by both observational tests is shaded. 



V. DISCUSSION 

We have performed a phase-space analysis on the dou- 
ble exponential potential model, its properties can be un- 
derstood in terms of a single exponential potential model 
with a varying coefficient A. The analysis shows that un- 
like the single exponential potential model which has two 
late times attractor depending on A, the double exponen- 
tial potential model has only one global attractor which 
leads to a matter dominated universe at late times. How- 
ever, it is always possible for the model to be dominated 
by the scalar field at present and dominated by matter 
in the future. The problem with this model is we need 
to fine-tune the parameters in order not to have reached 
the attractor at the present epoch, thereby circumvent- 
ing the objections of Binetruy [Q and de. la. Macorra 
|[13fl . However, all other quintessence models appear to 
also require a degree of fine-tuning 

We studied the observation constraints for both the 
single and double exponential potential models using a 
more updated type IA supernovae data and the gravita- 
tional lensing statistics. The results show that there are 
regions in the parameter spaces for which the models are 
consistent with the observations at the same time giving 
appropriate values for the scalar field energy density and 
the age of the universe at present. 
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